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ABSTRACT 


This  paper  deals  with  a  new  extension  of  a  weighted  residual  method 
of  analysis  called  the  extended  field  method.  The  extended  field  method 
Is  applied  for  the  first  time  to  the  problem  of  the  transient  vibration 
of  a  uniformly  thin  elastic  plate.  Numerical  results  have  been  obtained 
which  validate  the  analysis  procedure  and  show  better  solution  convergence 
than  Is  obtainable  by  standard  methods  of  analysis  for  the  same  number  of 
degrees  of  freedom.  Further  studies  are  necessary  to  obtain  still  better 
convergence. 


NOMENCLATURE 


a 

b 

c,c 

d 

e,e 

hj,n,k 


k,J.n.k 


k,m,n,p,q 

P(x,y;t) 

p 

Km,n 


extended  plate  length  In  x-coordlnate  direction 
extended  plate  length  In  y-coordlnate  direction 
extended  plate  offsets  In  x-coordlnate  direc¬ 
tion 

plate  thickness 

extended  plate  offsets  In  y-coordinate  direc¬ 
tion 

n,kth  Fourier  series  coefficient  of  the  extended 
plate  side  j  arbitrary  lateral  displacements, 
where  j  ■  1,2, 3,4 

n,k^  Fourier  series  coefficient  of  the  extended 
plate  side  j  arlbtrary  normal  slope  where 
j-  1,2, 3,4 
sunvnatlon  Indices 

pressure  loading  over  the  plate  area 

double  Fourier  series  coefficient  of  the  pressure 

loading 

farm — ->! 


spatially  uniform  pressure  amplitude 
time 

time  duration  of  pressure  pulse  loading 

plate  lateral  deflection  function 
particular  solution  component  of  w(x,y;t) 

plate  Cartesian  coordinates 

plate  stiffness  factor,  Ed"*/(12(l-u^) ) 

Young's  modulus 

maximum  values  for  the  indices 
1  if  i  is  odd;  0  otherwise 

constants  of  integration  for  complementary 
solution  to  the  non homogeneous  equation  (10) 
plate  lateral  deflection  amplitude;  the  sub¬ 
script  j  refers  to  a  specific  set  of  boundary 
conditions,  and  the  subscript  i  ranges  to  N 

(D/pd^Cfimr/a)^  +  (nir/b)2] 

(16Po0'm8n)/(^2'nnpd) 

Poisson's  ratio 
mass  density 

kir/T,  where  k  -  1 ,2 . N 


INTRODUCTION 


The  purpose  of  this  paper  is  to  present  an  extension  of  an  approxi¬ 
mate  method  of  analysis  called  the  extended  field  method  (XFM).  In 
brief,  the  XFM  employs  finite  series  which  are  term-by-term  solutions  to 
the  governing  differential  equation(s).  Then  the  device  of  an  offset 
auxiliary  boundary,  which  is  the  boundary  of  the  "extended  field",  is 
coupled  with  the  Galerkin  error  minimization  technique  in  order  to  ap¬ 
proximate  the  given  boundary  conditions  [1].  The  necessity  for 
term-by-term  series  solutions  to  the  governing  differential  equation 
limits  the  XFM  to  simple  geometries.  Despite  this  lack  of  versatility, 
the  XFM  is  noteworthy  because  of  its  demonstrated  superior  solution 
convergence. 

In  its  previous  development  the  XFM  has  been  applied  to  the  problem 
of  uniformly  thin,  elastic,  polygonal  plates  undergoing  forced  harmonic 
vibration,  and  small  combinations  of  plates  and  open  section  beams  un¬ 
dergoing  forced  harmonic  vibration.  For  the  purposes  of  discussing  the 
previous  displacement  amplitude  numerical  results,  single  plate  analyses 
are  divided  Into  two  groups.  The  first  group  consists  of  those  analyses 
where  the  polygonal  plates  possess  one  or  more  free  edges.  This  is  a 
difficult  group  for  the  XFM  because  the  free  edge  boundary  conditions 
require  three  differentiations  of  the  XFM  finite  series  solutions,  and 
there  Is  a  consequent  lessening  of  the  convergence  rate.  For  this  group 
of  plates,  as  many  as  ISO  degrees  of  freedom  may  be  necessary  to  achieve 
a  completely  convincing  four  significant  digits  of  convergence.  Even 
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so,  this  is  more  in  the  way  of  solution  convergence  than  can  be  obtained 
with  many  more  degrees  of  freedom  using  the  more  popular  methods  of 
analysis  such  as  the  finite  element,  finite  difference,  and  Rayleigh- 
Ritz  methods  [2]. 

The  second  group  of  analyses  concerns  polygonal  plates  without  free 
edge  boundary  conditions.  For  this  group  of  analyses,  the  XFM  generally 
attains  at  least  six  significant  figures  of  convergence  with  as  little 
as  one  hundred  degrees  of  freedom.  In  particularly  favorable  analyses, 
seven  significant  figures  of  convergence  can  be  attained  with  that  num¬ 
ber  of  degrees  of  freedom.  Thus  for  the  problem  of  linear,  forced  har¬ 
monically  vibrating  plates,  the  XFM  can  provide  benchmark  solutions. 
However,  forced  vibration  problems  are  a  relatively  simple  set  of  pro¬ 
blems. 

This  paper  presents  XFM  numerical  results  for  a  somewhat  more  com¬ 
plicated  problem,  that  of  the  linear  response  of  a  uniformly  thin  plate 
subjected  to  an  arbitrary  lateral  pressure  pulse.  In  this  problem  the 
time  variable,  t,  cannot  be  deleted  from  the  analysis,  and  initial  con¬ 
ditions  as  well  as  boundary  conditions  must  be  approximated.  Neverthe¬ 
less,  just  as  there  is  no  need  for  any  sort  of  geometric  grid  in  order 
to  describe  spatial  variations,  the  XFM  does  not  use  any  sort  of 
step-by-step  integration  process  to  evaluate  the  time  varying  response. 
Furthermore,  there  is  no  difference  in  style  between  this  XFM  analysis 
and  the  simpler  forced  harmonic  analysis  which  forms  a  portion  of  the 
present  analysis.  The  succeeding  sections  will  present  the  analysis 
procedure,  the  numerical  results,  and  the  conclusions  that  may  be  drawn. 


THE  ANALYSIS  PROCEDURE  IN  BRIEF 


In  Cartesian  coordinates  the  governing  differential  equation  for  a 
thin,  uniform,  isotropic,  homogeneous,  linearly  elastic  plate  subjected 
to  an  arbitrary  time  varying  pressure  loading  is 

DP4  w(x,y;t)+{pd)w>tt(x,y;t)  =  p(x,y;t)  (1) 

where  a  comma  indicates  partial  differentiation.  The  solution  to  this 
equation  is  obtained  in  several  parts  without  regard  to  the  shape  of  the 
actual  plate  boundary  or  the  nature  of  the  actual  boundary  supports. 
First  consider  the  homogeneous  form  of  Eq.  (1).  The  homogeneous  form 
can  be  solved  by  first  separating  the  spatial  and  temporal  variables  so 
that 

w(x,y;t)  *  W(x,y)sinwt  (2) 

Then  the  homogeneous  equation  becomes 

pfyx.y)  -  (pU2d/D)W(x,y)  *  0  (3) 

This  homogeneous  equation  is  the  same  as  the  homogeneous  equation  asso¬ 
ciated  with  forced  harmonic  vibration.  Thus  the  same  four  Levy  series 
solutions  used  In  the  harmonic  case  [1]  may  also  be  used  here.  However 
there  Is  the  difference  that  In  this  nonharmonic  case,  the  frequency  of 
forced  vibration,  u,  has  no  meaning.  In  order  to  determine  appropriate 
values  for  w,  define  a  time  period  of  Interest  (0,T)  over  which  the 


solution  will  be  sought.  In  engineering  terms,  this  time  period  is  sim¬ 
ply  defined.  It  would  usually  be  a  multiple  of  the  plate's  first  na¬ 
tural  period  or  a  multiple  of  some  characteristic  time  interval  associa¬ 
ted  with  the  applied  loading.  Therefore,  with  the  pulse  beginning  at 
t=  0,  write 

_  kl>  |.  _  l  o  1/  /  m  \ 


“  =  wk  "  T  * 


k  =  1,2 . K 


and  there  will  be  four  Levy  series  solutions  for  each  value  of  k.  Ig¬ 
noring  the  special  case  where  the  boundary  conditions  involve  adjacent 
free  edges,  which  is  dealt  with  in  Ref.  2,  the  solution  to  Eq.  (3)  is 

4  K 

w(x,y,t)  =  I  z  W.  t(x,y)sino).  t  (5) 

j=l  k=l  J,,<  K 


where 


W,>k(x,y)  ■  nriJ„,k{hl,n,ktsn,k  cos(h)sn,ka  sf"h  Vk<a-*> 

-rn,k  cosh  rn,ka  5f"0')sI,,k(a-»)3 
*  kl,n,k£s1"<h>  sn,ka  sinh  rn,k<*-*> 

-  sinh  rn  ka  sin(h)sn  k(a-x)]}sin 

Jn,k=  £sn,k  sinh  r"  tacos(h)s«  i  cosh  r„  ^  sin(h)s„  ,,a] 


n,k  n,k 


rn,k  =  Knn/b: >  +  Wk(pd/D)^ 
sn,k  *  [{"’■/t.2)  -  ^(pd/D)*5]** 
-  [™k(pd/D)!i  -  (nn/b)2]’1 
r>k  =  {b/*)(pukd/D)'s 

s1"<h>sn,kx  ■  s<"5„,kx 
s(nhsnkx 
cos(l.)sn>kX  »  cos  s„  kx 
cosh  s„  kx 


A 
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for  n  <  n. 


if  n  < 

A 

if  n  >  nk 

A* 

if  n  <  nk 

A 
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As  a  brief  explanation,  for  example,  the  Levy  series  k  can  be  thought 

of  as  representing  the  solution  for  a  rectangular  plate  of  dimensions  a 
by  b  with  simple  supports  at  sides  y*0,b  and  x=a,  but  arbitrary  la¬ 
teral  displacements  and  normal  slopes  at  side  one  (x*0)  that  are  de¬ 
scribed  by  the  respective  equations 

N  K 

H(°,y)  -  l  I  h1  n  k  sin  ^ 
n-1  k*l  l,n,K  D 
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(7) 


W  JO.y) 


N  K 
r  z 
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1  .n,k 


This  rectangular  plate  is  called  the  extended  plate.  Later  in  the  anal¬ 
ysis,  the  rectangular  extended , plate  boundaries  will  enclose  the  actual 
plate  boundaries.  The  third  Levy  series,  n  ^ ,  which  corresponds  to 

arbitrary  boundary  conditions  on  the  extended  plate  edge  x=a,  is  like 
the  first  series,  but  with  x  replacing  a-x  and  -kj  n  k  replacing 

♦ki  jj.  The  second  series,  corresponding  to  arbitrary  conditions  at 

y*0  can  be  obtained  from  the  first  by  interchanging  a,x  and  b,y.  The 
fourth  solution  can  be  obtained  from  the  second  as  the  third  is  obtained 
from  the  first. 


To  obtain  the  particular  solution  for  the  non-homogeneous  equation 
(1)  write  the  forcing  pressure  in  double  Fourier  series  form  as 

M  M 

p(x,y;t)  -  z  z  pm  n(t)sin(mirx/a)sin(mry/b)  (8) 

m=l  n=l  m,n 


where  the  index  M  is  large  enough  to  ensure  satisfactory  convergence. 
Then  the  particular  solution  may  be  written  in  Navier  series  form  as 

M  M 

w5(x,y,t)  =  Z  Z  Zm  n(t)sin(miTX/a)sin(nity/b)  (9) 

m=l  n=l  * 


Substitution  of  Eq.  (9)  into  Eq.  (1),  and  noting  the  linear  independence 
of  the  sine  functions,  leads  to 


(t)  + 


Wtl  "  Wt,/Pd 


(10) 


wh6r6  a  p  op 

a  =  (D/pd)](imt/a)  +  (n-rr/b)  ]  . 

Ill  )  fl 
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The  solution  of  these  M  elementary  ordinary  differential  equations  com 
pletes  the  solution  for  Wg(x,y,t),  and  hence  for  w(x,y,t)  which,  again, 

is  the  right  hand  side  of  Eq.  (5)  plus  Wg(x,y;t). 


In  order  to  illustrate  the  above,  consider  a  fully  clamped  rectan¬ 


gular  plate  of  dimensions  a-c-c  by  b-e-e  where  c,e,c,e  are  the  offsets 
between  the  actual  plate  rectangular  boundary  and  the  extended  rectan¬ 
gular  plate  boundary.  Let  the  applied  loading,  which  is  extended  to  the 
boundaries  of  the  extended  plate  as  well,  be 

0  <  t  <  t_ 


p(x»y»t) 


( 


pQ  s1n{*t/t0) 


-  o 
otherwise 


(ID 


This  choice  of  plate  and  loading  for  an  example  problem  is  neither  fa¬ 
vorable  or  unfavorable  for  the  analysis  procedure.  The  four  clamped 
edges  require  a  larger  number  of  degrees  of  freedom  than  is  necessary 
In  many  cases,  and  the  double  Fourier  series  representation  of  p(x,y;t) 
is  slow  to  converge.  On  the  other  hand,  the  absence  of  one  or  more 
free  edges  Improves  the  convergence  rate  [2].  Substitution  of  Eq.  (11) 
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into  Eq.  (8)  yields 


where  O' 
m 

solution 


pm,n{t)  =  (16PoflrmV,r2mn)sin^t/to) 
is  plus  one  if  m  is  an  odd  integer,  and  zero  if  even, 
to  Eq.  (10)  is  easily  found  to  be  such  that,  for  0  <  t 


w5(x,y,t) 


M  M 
z  z 
m=l  n=l 


wm,n' 


m,n 


(12) 
Then  the 


m,n 


a  -(u/t  ) 
m,n  x  '  o7 


j  sin(7rt/t0)}sin(m7rx/a)sin(n7ry/b) 


(13) 
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where  &  „  =  1 6p^GT  G^/tt  mnpd  and  Qm  „  and  R  „  are  constants  of  integra- 
m,n  ro  m  n  Tn,n  m,n  3 

tion.  The  time  period  after  t =  t  can  be  dealt  with  by  standard  super¬ 
position  procedure  involving  two  continuous  sinusoidal  loadings  with  a 
half  period  phase  shift. 


Now  it  is  necessary  to  satisfy  the  initial  conditions  and  the  ac¬ 
tual  plate  boundary  conditions.  The  initial  conditions  will  be  satis¬ 
fied  first.  Let  these  conditions  be  those  of  zero  initial  displacement 
and  zero  initial  velocity.  Zero  initial  displacement  quickly  leads  to 
Rm  n=0.  To  satisfy  zero  initial  velocity,  K  is  set  equal  to  N,  and 

Qm  n  is  split  into  two  parts.  Q-j  m  n  is  used  to  cancel  the  w&  portion 

of  the  initial  velocity  equation,  and  m  n  cancels  the  through  W4 

portion.  In  other  words,  due  to  the  linear  independence  of  the  sine 
functions 


"n3_  Ja 


M  ,m,n 


m^rr 


mjn  o 


am,n  “  (ff/to>2 


m,n  =  1 ,2,. . .  ,M 


(14) 


and  Q2  m  n  satisfies  the  equation 
N  N 

1  1  Q?  m  m  “m  n  sin(mnx/a )sin(niry/b) 

m=l  n=l  m,n 


k=1 


V-tu 


l,k 


W4,k^  = 


(15) 

(16) 


where  Q2  m  n  is  zero  for  m,n  greater  than  N.  (Note  M»N.)  Eq.  (15)  is 

written  as  an  equality,  but  it  is  actually  only  an  approximation.  The 
error  of  this  approximation  may  be  minimized  by  use  of  the  Galerkin 
technique  with  weighting  functions  sin(pitx/a)  and  sin(qny/b)  on  the  in¬ 
tervals  (0,a)  and  (0,b).  For  the  sake  of  brevity,  the  algebraic  solu¬ 
tion  for  Qo  „  „  will  be  omitted  here. 
y  ,m,n 

The  satisfaction  of  the  eight  boundary  conditions  is  even  more 
Straight  forward.  For  example,  the  boundary  condition  w(c,y;t)  =  0  is 
approximated  by  substitution  of  the  expression  obtained  above  for 
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w(x,y;t).  Then  the  error  of  that  approximation  to  the  boundary  condi¬ 
tion  w(c,y,t)  =  0  is  minimized  by  the  Galerkin  technique  along  the  length 
of  the  actual  edge  x  =  c  by  use  of  the  weighting  function  sin(piry/b)  on 
the  interval  (0,b)  and  minimized  again  over  the  time  interval  (Q,T)  with 
the  weighting  function  sin  (q-rrt/T ) ,  where  p,q  =  1,2,...,N.  The  result  is 

2  2 
N  linear,  non-homogeneous  algebraic  equations  in  the  8N  unknown  ex¬ 
tended  plate  edge  coefficients  h^  ^  ^  through  k^  N  The  other  seven 

o 

boundary  condition  equations  supply  the  remaining  7N  equations.  Ref.  3 
details  the  above  derivation  and  lists  all  necessary  results,  specifi¬ 
cally,  the  unabridged  forms  for  the  eight  boundary  condition  equations 
whose  solution  explicitly  defines  w(x,y;t). 


NUMERICAL  RESULTS 


The  fully  clamped  rectangular  plate  subjected  to  the  half-period 
sine  pulse  discussed  above  was  programmed  for  the  UNIVAC  1140  digital 
computer.  The  selected  input  for  the  plate  geometry  and  material  pro¬ 
perties  was  a  =  20  inches,  b  =  16  inches,  d  =  1  inch,  c  =  e  =  c  =  e=0.0, 

2  4 

E=  10,500,000  psi,  p  =  0.30,  and  p  =  0.000254  lb-sec  /in  .  The  selected 
parameters  for  the  pressure  loading  were  a  uniform  spatial  amplitude  of 
pQ=l  psi,  and  a  half  sine  pulse  duration  of  tQ=  0.010873  sec.  The 

latter  number,  besides  having  too  many  digits  for  an  essentially  arbi¬ 
trary  choice,  was  based  on  an  originally  mistaken  estimate  of  the 
plate's  first  natural  period.  A  correct  estimate  is  0.00087  secs.  [4]. 
Thus  the  selected  pulse  duration  turned  out  to  be  approximately  twelve 
and  one-half  times  the  plate's  first  natural  period.  Thus  it  was  rea¬ 
sonable  to  set  the  time  period  of  interest,  T=tQ.  In  order  to  well  de¬ 
fine  the  response  in  the  interval  (0,T),  it  was  decided  to  calculate  the 
displacement  response  at  101  equally  spaced  points,  including  the  end 
points. 

The  primary  purpose  for  gathering  the  numerical  results  presented 
here  was  to  investigate  the  degree  of  solution  convergence  obtainable  by 
the  XFM  analysis,  and  thus  learn  if  the  same  high  degree  of  convergence 
that  characterized  forced  harmonic  vibration  XFM  solutions  could  also 
be  obtained  for  pulse  loadings.  In  order  to  study  the  XFM  solution 
convergence,  three  indices  were  varied.  The  first  of  these  was  N,  the 
maximum  number  of  terms  in  each  solution  series.  The  importance  of  N 
is  that  the  number  of  degrees  of  freedom  that  describe  the  plate  dis- 
2 

placements  is  8N  .  This  index,  as  expected,  was  the  critical  index. 

The  second  index  was  N,  another  maximum  summation  index  originally  equal 
to  N,  but  by  rearranging  certain  orders  of  summation,  could  be  freed 
from  being  tied  to  the  small  value  that  N  must  normally  be.  The  results 

in  Table  1  show  that  a  value  of  N*25  is  quite  sufficient.  The  third 
index  that  was  varied  is  M,  which  determines  the  number  of  series  terms 
used  to  describe  the  input  pressure.  As  Table  1  suggests,  M,  which  is 
an  odd  number  because  even  indexed  input  terms  are  zero,  is  less 
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influential  than  N.  In  the  two  instances  checked,  only  the  sixth  sig¬ 
nificant  digit  was  affected  when  M  was  increased  from  49  to  99,  and  in 
the  one  instance  checked,  only  the  seventh  significant  digit  was  affec¬ 
ted  when  M  was  increased  from  99  to  149.  Thus  only  N  needs  further  con¬ 
sideration  when  discussing  convergence  at  this  time. 

Before  studying  the  effect  of  N  on  convergence  it  is  worth  noting 
that  the  digital  computer  program  that  supplied  the  results  presented 
in  Table  1  was  verified  by  comparison  to  a  finite  element,  direct  step¬ 
wise  integration  solution.  Specifically,  an  XFM  solution  at  N  =  8, 

N  =  49,  M=99  for  the  same  plate,  but  for  a  change  in  thickness  to  0.1 
inch,  was  within  one  half  percent  of  the  corresponding  NASTRAN  solution 
using  eighty  QUAD2  elements  and  364  degrees  of  freedom  over  one  quarter 
of  the  doubly  symmetric,  actual  plate.  This  is  a  good  place  to 
mention  that  although  the  XFM  solution  could  also  have  been  greatly  sim¬ 
plified  by  use  of  the  double  symmetry  of  the  plate  geometry  and  loading, 
no  such  simplification  was  made  in  order  to  present  results  that  equally 
well  characterize  a  loading  with  any  spatial  distribution. 

The  digital  listing  of  the  time  response  curve  for  a  particular 
geometric  point  on  the  plate  represents  a  lot  of  different  numbers.  Re¬ 
sponse  curves  for  different  maximum  indices  could  have  been  compared  on 
the  basis,  for  example,  of  an  averaged  sum  of  the  absolute  values  of  the 
displacement  differences,  or  some  similar  scheme  that  compares  the  en¬ 
tire  curves  on  a  weighted  or  unweighted  basis.  However,  from  an  engi¬ 
neering  point  of  view,  the  maximum  displacement  is  usually  the  most  in¬ 
teresting  displacement  measurement.  Of  course,  the  true  maximum  lies 
between  data  points.  Table  1  presents  the  values  at  that  data  point 
(time  instant)  closest  to  where  the  true  maximum  occurs.  In  all  cases 
studied  it  was  the  same  time  instant,  number  56,  which  is  equivalent  to 
an  elapsed  time  of  0.00598015  sec.  Returning  now  to  the  convergence 
pattern  dependent  upon  N.  It  suffices  to  look  at  the  data  for  which 

M=99  and  N=49.  The  plate  center  displacements  for  N=3,4,  and  5  have 
only  converged  to  the  third  significant  digit.  This  degree  of  conver¬ 
gence  is  typical  of  the  convergence  at  other  time  points,  with  very 
little  deviation  in  the  convergence  rate  from  one  time  point  to  another. 
This  situation  is  to  be  expected  since  the  solution  error  is  essentially 
the  same  at  all  instants  of  time.  That  is,  the  error  does  not  increase 
with  time  as  it  does  with  a  step-by-step  integration  scheme.  The  same 
background  comments  apply  for  the  values  of  N  equal  to  six  and  seven. 
Here  it  is  almost  possible  to  claim  a  fourth  digit  of  convergence  be¬ 
cause  the  Cauchy  differences  are  decreasing  rapidly.  Specifically,  the 
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Cauchy  differences  multiplied  by  10  for  the  pairs  of  values  of  N  equal 
to  (5,4),  (6,5),  and  (7,6)  are  67,  50,  and  5  respectively.  However,  it 
is  not  possible  to  guarantee  that  the  above  trend  will  continue.  There¬ 
fore  only  convergence  to  a  third  digit  will  be  claimed,  or  in  other 
terms,  the  maximum  center  displacement  appears  to  be  0.1304  ±  .0004  x 

10'3  inches. 


CONCLUSIONS  AND  FUTURE  WORK 


From  the  comparison  to  the  NASTRAN  result  it  is  clear  that  the  XFM 
does  provide  correct  solutions  for  the  pulse  loaded  thin  plate  problem. 
The  three  significant  digit  convergence  at  two  hundred  degrees  of  free¬ 
dom  is  not  spectacular,  but  it  is  still  distinctly  better  than  can  be 
obtained  from  a  standard  finite  element,  finite  difference,  or  Rayleigh- 
Ritz  analysis.  It  is  also  clear  that  solutions  for  increased  values  of 
N  are  required  to  more  fully  explore  the  XFM  convergence  rate.  Other 
parametric  studies  are  also  needed.  Variations  on  the  relations  between 
the  pulse  duration,  the  first  natural  period,  and  the  period  of  interest 
need  investigation.  Reasonable  variations  on  the  offset  distances  and 
pulse  shape  are  not  expected  to  have  any  effect  on  solution  convergence 
but  this  also  should  be  verified. 
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